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Abstract 

The Universal Extra Dimensions model has been implemented in the Pythia generator from 
version 6.4.18 onwards, in its minimal formulation with one TeV -1 sized extra dimension. The 
h— j | additional possibility of gravity-mediated decays, through a variable number of eV _1 sized extra 

\^ ■ dimensions into which only gravity extends, is also available. The implementation covers the 

fSJ | lowest-lying Kaluza-Klein (KK) excitations of Standard Model particles, except for the excita- 

tions of the Higgs fields, with the mass spectrum calculated at one loop. 2^2 tree-level pro- 
' ' duction cross sections and KK number conserving 2-body decays are included. Mixing between 

iso-doublet and -singlet KK excitations is neglected thus far, and is expected to be negligible for 
■ ' all but the top sector. 

1 Introduction 

| In the Universal Extra Dimensions (UED) model, first formulated in [1], all Standard Model (SM) 

fields are allowed to propagate into 5 TeV -1 sized extra dimensions. In its minimal formulation, 
the SM lives in 4 + 5(= 1) space-time dimensions. This model can be considered as an effective 
theory, valid below some cutoff scale A > l/R, where R is the compactincation length of the extra 
dimension. To avoid fine-tuning of the parameters in the Higgs sector, l/R should also not be much 
higher than the electroweak scale. Such models can be shown to be consistent with all current low- 
energy and collider constraints [1-3]. 

If the UED space is further embedded into a larger space with N extra dimensions into which 
^ ■ only gravity spreads [4], gravity mediated decays also become possible. In this case the (4 + AO- 

dimensional Planck scale Md should not be more than one or two orders of magnitude above l/R [5]. 

Phenomenologically, UED models exhibit several interesting properties, often similar to those of 
supersymmetric (SUSY) models. Every SM field has a Kaluza-Klein (KK) partner (a whole tower 
of them in fact, but here we consider only the lowest-lying excitations), which carries a conserved 
quantum number, KK parity. This conserved parity, tracing its origin to extra-dimensional momentum 
conservation, renders the lightest KK particle (LKP) stable. Heavier KK modes cascade decay to 
the LKP by emitting relatively soft SM particles. The LKP escapes detection, generally resulting 
in missing energy signals. If the model is extended to include (N + 4)-dimensional gravity, then 
the LKP may decay further to its SM partner plus a low-mass graviton excitation, in which case 
a phenomenology more similar to that of SUSY models with gauge mediated SUSY breaking may 
arise [6,7]. 

In the Pythia [8] implementation, the adjustable parameters for the default UED model (i.e., with 
no modification of the gravitational sector) are l/R, A, and the number of quark flavours in the KK 
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excitation spectrum, with default values of 1/R = 1 TeV, A = 20 TeV, and = 5, respectively, 
while 5 is fixed to the value of 1. Optionally, one may choose to adjust the value of A x R rather than 
that of A. For the variant of the model with (4 + AQ-dimensional gravity, the number of eV _1 -sized 
extra dimensions, N = 2, 4, or 6, and the scale Md also enter as adjustable parameters, with default 
values of A^ = 6 and Mq = 5 TeV, respectively. 

In Section [2 we give a brief description of the spectroscopy of the minimal UED model and in 
Section [3] extend it to include gravity mediated decays. In Section 01 the UED implementation in the 
Pythia generator is detailed. We conclude and give an outlook in Section [5] 



2 One-dimensional UED spectroscopy 

In one-dimensional UED, each SM particle has n = 1, 2, 3, ... KK excitations, of squared mass 



m, 



m| M + n 2 /R 2 



(1) 



with the n = state corresponding to the SM particle. So far, only the n = 1 excitations have been 
incorporated in Pythia. The compactification of the extra dimensions is assumed to conserve the 
extra-dimensional momentum components at tree level, resulting in conservation of KK number in 
the effective 4-dimensional model. Hence KK particles are always produced in pairs and there are no 
vertices involving only one non-zero KK mode. 

At tree level, the mass spectrum at each KK level is highly degenerate for all SM particles with 
msM C cf. eq. (Q]). However, first order radiative corrections, calculated in [9] and implemented 
in Pythia using the code from [10], lift this degeneracy by about 20% for strongly interacting par- 
ticles, the heaviest being the excited gluon, and by less than 10% for leptons and EW bosons, the 
lightest being the excited photon (the LKP). With these mass splittings, the SM quark and gluon KK 
excitations cascade decay down to the LKP, which is stable (unless (N + 4) -dimensional gravity is 
switched on, see below). An example of a nondegenerate KK particle mass spectrum obtained with 
Pythia for representative values of the free parameters is shown in Figure Q] These results agree 
with [9]. 
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Figure 1: In the first figure is shown the KK particle mass spectrum of the first level KK states including radiative 
corrections [9], for 1/R = 500 GeV, AR — 20 and a„ = 0.118. The star (*) denotes a KK particle and the notations D 
and 5 indicate respectively the doublet and singlet KK fermions. In the second figure are shown the dominant (solid) and 
rare (dotted) decays of KK particles. The particles not denoted by a star (*) are SM. 



2 



In general, mixing can occur between the doublet and singlet KK states. This effect is strongly 
suppressed for the light-flavor KK states, but can be phenomenologically relevant for the KK top 
quark [1]. Nonetheless, since the excited top pair production is small compared to the sum of all 
KK processes at the LHC [11], neglecting the mixing should still be a reasonable approximation for 
observables that are not explicitly sensitive to the top flavor. In the current Pythia implementation, 
the effects of doublet-singlet mixing in the KK sector are thus neglected. 

3 Gravity mediated UED decay widths 

If the (4 + l)-dimensional UED space (brane) is embedded into a larger space of (4 + N) dimensions 
(bulk), where iV counts the number of eV _1 sized extra dimensions into which only gravity propagates 
(with one of the N being parallel to the UED dimension), then gravity mediated decays also become 
possible [4]. The phenomenology of these decays then also depends on the (4 + iV)-dimensional 
Planck scale Mrj, introducing two additional free parameters in the model, N and Mr> 

The graviton field appears as a massless particle with a tower of excited modes whose masses dif- 
fer by order of eV. Most importantly, the graviton modes extending in the UED direction couple to KK 
particles. Hence KK particles may decay directly to SM particles by emitting such low-mass graviton 
excitations. While the coupling to each such mode is incredibly small, the modes are sufficiently 
densely spaced that, after summation (or, in the continuum approximation used for practical calcula- 
tions, integration) over them, total transition rates relevant for collider phenomenology may occur. In 
general, the gravity-mediated decays will compete with the non-gravitational (mass-splitting) modes. 
If 

r(mass splitting) > T (gravity mediated) 

then the gluon and quark excitations will cascade down to the excited photon 7* (LKP), which then 
will decay via 7* — > 7 + G^*\ More details can be found in [12]. This decay is the only gravity- 
mediated mode that appears by default in the Pythia implementation when (N + 4) -dimensional 
gravity is switched on. The branching ratio for 7* — > 7 + G^ is then 100%. 

The remaining gravity mediated decays (of all other KK particles) are foreseen to be included in 
a future version and are currently available as a standalone add-on routine, pygrav . f , which can be 
downloaded from [14], with width expressions from [5, 13, 15]. 

The graviton mass in these decays is obtained by integrating the differential width taken fromQ [5, 
13]. The formulae for the mass splitting decay widths in PYTHIA were extracted from the code [10] 
and their dependence on 1/Ris illustrated in Figure |2j which is in agreement with [12]. 

4 Pythia implementation 

The Pythia implementation of the UED particle spectrum, production processes, and decay modes 
is summarized here, after which we give a brief overview of the relevant user switches, parameters, 
and subroutines controlling the code. 

4.1 Pythia UED particle spectrum, production processes and decay modes 

In Pythia, in order to avoid confusion between chiral and weak eigenstates, the UED KK states are 
labeled following a slightly different convention than that so far adopted by the PDG. In the PDG, the 

'This corrects a numerical problem in the implementation of [1 1]. 
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Pythia 6.4.20 UED - 1st level KK fermions 
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Figure 2: On the left-hand figure, the mass splitting decay widths are shown versus 1/72 for the first level KK excitations 
of vector bosons for the following decays: g* — > ddj), W** — > e^D an d Z* — > e + ej} _ . On the right-hand figure, the 
mass splitting decay widths are shown for the first level KK excitations of fermions : qjj, q s (for up-type quarks; the width 
for down-type quarks is four times smaller) and These plots were produced using PYTHIA 6.4.20 with AR = 20. 



KK particles are labeled such that for the fermions, the subscripts L and R denote that the fermion is 
respectively a doublet or a singlet under SU(2)w (see Table Q] for first level excitations). However, 
these doublet and singlet UED excitations are ordinary Dirac fermions, which both have left- and 
right- handed chiral spinor components. 

To avoid any confusion with helicity states, the notations D (for doublet) and S (for singlet) are 
used here instead of the PDG L and R. Furthermore, the superscript W is replaced by a star (*) to 
correspond to the usual Pythia notation for extra-dimensional excitations. The relationship between 
the Pythia and PDG particle names is given in Table Q] We emphasize that these are only notational 
differences, of a purely cosmetic nature. We use the same numbers (particle codes) as the PDG, and 
our scheme is therefore fully compatible with theirs. 

The UED states can be produced through nine new production processes (ISUB=311 to 319), 
listed in Table |2] These employ tree-level differential cross section expressions [5, 11] and their 
dependence on 1/R is illustrated in the two plots of Figure [3] (for default choices of all other Pythia 
parameters, such as the strong coupling, parton distributions, and renormalization and factorization 
scales). These can be compared to those in [5, 12]. UED processes generated using an external 
generator can of course also be interfaced, using the existing LHA [16] and LHEF [17, 18] interfaces, 
in which case Pythia will handle the subsequent decays, radiation, and fragmentation. 

Finally, the mass-splitting decay modes as well as their typical branching ratios and widths are 
given in Table [3] Again, externally calculated branching ratios can also be interfaced, if so desired, 
using the existing interface to SLHA decay tables [19]. As mentioned previously, in this version of 
Pythia, the excited states of the Higgs are not implemented, and the only available gravity mediated 
decay is that of the excited photon, 7* — ► 7 + G^*\ 
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Table 1 : PDG and Pythia notations and codes for the first level KK excitations. 



ISUB 


Production process 


Note 


311 


g + g -> g* + g* 




312 


g + q^g* + qo5 g*+qg 




313 


qi + qj ^qDi + qiV q^ + qsj 


all i, j 


314 


g + g^qo + qo; qs + ^s 




315 


q + q^qD + qD; qs + qs 




316 


cu + qj -> qfji + qsj 


« / J 


317 


qi + qj ^qfji + qDj; qsi + qsj 




318 


cu + qj -» qDi + 


all i, j 


319 


Oi + Qi -> q£)i + qioi 


all i, j 



Table 2: UED production processes and their associated PYTHIA ISUB process number. 

4.2 Pythia user switches and subroutines for UED 

The UED parameters which can be modified by the user are: 

• the compactification scale or curvature of the extra dimension, 1/R, 

• the cutoff scale of the theory, A (or, alternatively, A x R), 

• the number of quark flavors, 

• whether to use the extension to (4 + N) -dimensional gravity and hence allow LKP decay by 
graviton emission, and if so, 

• the number of large extra dimensions where only the graviton propagates, N, 

• and the (4 + AO-dimensional Planck scale, Md, 

The Higgs boson mass is also a free parameter in the UED theory but it is set through the usual 
Pythia pma s (25, 1 ) parameter. In the code, the UED switches and parameters are stored in the 
new common block: 
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Figure 3: Cross sections for the production of two stable KK final states in proton-proton collisions at E cm = 14 
TeV (LHC), generated using PYTHIA 6.4.20. On the left-hand figure are shown the cross sections versus 1/R for quasi 
degenerate KK particle masses ttikk — 1/R (user switch IUED(6)=0; see Section l4~2t . for the different production sources: 
gg, gq and qq, and for the sum of the three. On the right-hand figure are shown the cross sections for the production of 
two hard photons with missing transverse energy in the final state, where KK particle masses include radiative corrections 
(IUED(6)=1, the user switch default value), and with the following kinematic cuts : P^SPt > 200 GeV/c 2 and E™ iss > 
200 GeV, for N = 2 and 6. For all cases, AR = 20. 



COMMON/PYPUED/IUED (0:99) , RUED (0:99) 

IUED ( 1 ) = The main UED ON(=1)/OFF(=0) switch 
Default value = 

IUED ( 2 ) = On/Off switch for the extension to (N + 4)-dimensional gravity (switching it 
on enables gravity-mediated LKP decay): ON(=1)/OFF(=0) 
Default value = 

IUED ( 3 ) = The number of KK excitation quark flavors 
Default value = 5 

IUED ( 4 ) = N, the number of large extra dimensions where only the graviton propagates. 
Only used when IUED (2) =1. 
Default value = 6 (can be set to 2, 4 or 6) 

IUED (5) = Selects whether the code takes A (=0) or AR (=1) as input. See also 
RUED (2:3). 
Default value = 

IUED ( 6) = Selects whether the KK particle masses include radiative corrections (=1) or 
are nearly degenerate ttikk — 1/R (=0). 
Default value = 1 
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Decay mode 


Branching ratio 


Total Width (GeV) 
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1/3 




— * Qi + W* (except for q = b, t) 


2/3 




b + Z* 


100% 


5.21 x 10~ 2 


t* _> b + w* 


100% 


1.04 x 10 _1 




1 /fi 

1/U 


% 7^ v i n _1 


Z* -» i/ + i£ (and i/ + pft) (and P + l*,*) 


1/12 


1.92 x 10 -1 


g* -> all 


100% 


53.9 


g* — > q + q§ (and q§ + q) (down type) 


6.4% 




g* — > q + q§ (and qg + q) (up type except for q = t) 


6.0% 




g* — > q + qjp (and qj^ + q) (except for q = t) 


3.8% 





Table 3: UED decay modes, branching ratios and widths for 1/R = 500 GeV and Ai? = 20. Note that 
certain decays are kinematically suppressed for lower values of 1/R. This is the case when 1/R = 500 
GeV for t§ -» t + 7*, b^ -» t + W*, tf, -► t + Z* and g* -► t + t§ or D . The q* s -» Z* + g and 
Z* — > /J + / decays have been switched off due to the fact that they are suppressed by a factor sin 2 6\, 
the level 1 Weinberg angle which is of order 10~ 2 — 10 -3 , whereas the W* — > l* s + v decay is 
forbidden [20]. 



RUED ( 1 ) = 1/R, the curvature of the extra dimension 
Default value = 1000 GeV 

RUED ( 2 ) = M D , the (4 + AO-dimensional Planck scale. Only used when I UED ( 2 ) =1. 
Default value = 5000 GeV 

RUED ( 3 ) = A, the cutoff scale. Used when IUED ( 5 ) =0. 
Default value = 20000 GeV 

RUED ( 4 ) = Ai?, the cutoff scale times the radius of the extra dimension. Used when 
IUED (5) =1. 
Default value = 20 



Four new subroutines and two new functions were added to handle UED-specific tasks: 

SUBROUTINE PYXDIN to initialize Universal Extra Dimensions 

SUBROUTINE PYUEDC to compute UED mass radiative corrections 

SUBROUTINE PYXUED to compute UED cross sections 

SUBROUTINE PYGRAM to generate the UED KK graviton mass spectrum 

FUNCTION PYGRAW to compute UED partial widths to G* 

FUNCTION PYWDKK to compute UED differential widths to G* 
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In addition, several Pythia routines were modified to facilitate the UED implementation. These are 
SUBROUTINE PYGIVE now accepts input also for IUED and RUED 
SUBROUT INE PYINIT added call to PYXDIN to initialize UED 
SUBROUTINE PYMAXI small extension for UED overestimates 
SUBROUTINE PYPTFS small extension for showering KK gluons 
SUBROUTINE PYRAND extended to choose flavors in UED processes 
SUBROUTINE PYRESD added call to PYGRAM to choose graviton mass from 

continuous spectrum in UED decays to gravitons 
SUBROUTINE PYSCAT extended to include UED processes 
SUBROUTINE PYSIGH small extension to call PYXUED for UED 
SUBROUTINE PYWIDT extended to compute KK decay widths 

5 Conclusion and outlook 

The minimal UED (mUED) model with one extra dimension has been implemented in the Pythia 
generator from version 6.4. 18 onwards. The additional possibility of gravity mediated decays has also 
been included. The model uses 1-loop corrected mass formulae and tree-level expressions for cross 
section and decay width calculations. 

The main point of this work is to facilitate complete collider phenomenology studies of UED 
signatures, by combining the leading-order production and decay matrix elements discussed in the 
main body of thispaper with the more traditional components of the Pythia generator: sequential 
resonance decays q parton showers, hadronization, and modeling of the underlying event. Due to the 
typically large absolute mass scales and relatively small mass differences in the mUED model, addi- 
tional QCD jets from initial-state radiation can be an important source of combinatorial error when 
attempting to identify the jets emitted in decays of colored KK particles H The UED implementa- 
tion is compatible with both the old (Q 2 -ordered [23-25]) and new (p\ -ordered [26]) shower and 
underlying-event models available in Pythia 6.4. Hopefully, this will make it easier to evaluate not 
only the overall impact of the QCD collections, but also to gain some insight into their uncertainties. 
Large QCD uncertainties may sometimes be reduced by the application of matrix-element-to-parton- 
shower matching methods (see, e.g., the reviews in [22, 27-29]), but this was deemed beyond the 
scope of the present work. 

This work was started at the Les Houches Workshop in 2005 [30] first using the CalcHEP and 
CompHEP event generators [31] and since then ongoing work has been carried out in the ATLAS [32- 
34], International GDR [35], and MC4BSM [36] contexts. 

The next step will be to implement the model in C++ in Pythia 8 [37]. We furthermore intend to 
include the whole set of gravity mediated decay widths, but deemed this low priority at present, since 
only a small region of parameter space is concerned, where the mass splitting and gravity mediated 
widths are of the same order of magnitude. Likewise, for studies concentrating on the top sector, the 
doublet-singlet mixing effects in the KK top sector should be included. Finally, it would be interesting 



"As an alternative to the decays included in this implementation, externally generated decay tables for the UED particles 
can also be read in, e.g., using the SLHA format [19]. 

3 For illustration, see, e.g., the corresponding case for SUSY, studied in [21, 22]. 
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to include also the excitations of the SM Higgs fields as well as the effects of potentially resonant KK 
number violating interactions mediated by the 2nd level KK states [9,20, 38]. 
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